Probing the local temperature of a 2DEG microdomain with a quantum dot: 
measurement of electron-phonon interaction 
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We demonstrate local detection of the electron temperature in a two-dimensional microdomain 
using a quantum dot. Our method relies on the observation that a temperature bias across the dot 
changes the functional form of Coulomb-blockade peaks. We apply our results to the investigation 
of electron-energy relaxation at subkelvin temperatures, find that the energy flux from electrons into 
phonons is proportional to the fifth power of temperature and give a measurement of the coupling 
constant. 



o 

(N 

=5 



I 



a 



I 

o 
o 



> 

1^ 



1^ 

O 
O 

> 

X 



The role of temperature and the investigation of heat 
transport in solid-state nanoscale systems are currently 
under the spotlight [T]. Electronic thermometers [2J- 
H], refrigerators [5HH] and heat transistors [5] based on 
metallic and semiconductor nanostructures were recently 
demonstrated. Such devices are typically operated at 
subkelvin temperatures, where electron-phonon interac- 
tion is weak and, as a result, the relevant electronic tem- 
perature may significantly differ from that of the lattice. 
In this context, the ability to perform local measurements 
of the electronic temperature becomes highly desirable. 
While in metallic systems superconducting tunnel junc- 
tions are routinely employed to this end, for the case of 
two-dimensional electron gases (2DEGs) or semiconduc- 
tor nanowires it was shown that quantum dots (QDs) can 
serve the purpose |Tni |TT] , due to the fact that at thermal 
equilibrium the linewidth of zero-bias Coulomb blockade 
(CB) peaks in the weak-coupling regime is directly re- 
lated to the electronic temperature of the leads. 

In this Letter we propose a method for the detection 
of local temperature in a 2DEG based on the analysis 
of the CB-peak lineshape in the presence of a temper- 
ature bias, and apply it to the investigation of energy- 
relaxation mechanisms in a micrometer-sized electronic 
domain. The domain is defined electrostatically and con- 
nected to the surrounding 2DEG regions through a QD 
and three quantum point contacts (QPCs). As a known 
heating power is delivered to the domain, conductance 
across the QD is probed and the steady-state tempera- 
ture detected. The dependence of the measured electron 
temperature on the heating power is studied for differ- 
ent values of the QPC resistances: as the coupling to 
the surrounding 2DEG regions is reduced, we observe 
the crossover from a regime where excess heat is carried 
away by hot quasiparticles tunneling through the QPCs 
to one where power exchange with lattice phonons dom- 
inates. Our measurements show that the latter mech- 
anism follows the T^ power law expected [12] for the 
screened electron-acoustic phonon piezoelectric interac- 
tion and yielded a measurement of the coupling constant. 

Figure 1(a) shows a scanning electron micrograph 



(a) 



Vac (y^. 



(b) 0. 050 G (mS ) 5.3 25 
-0.40 , 




-1.0 



Vopc, (V) 



20 -10 10 20 
l(nA) 



FIG. 1. (color online), (a) False-color scanning electron 
micrograph of the device and scheme of the measurement 
setup. Schottky gates (light gray) define a IGjini'^ sized cen- 
tral electronic domain connected to outer regions of the 2DEG 
through a QD (bottom) and three QPCs (left, right, top). 
Crossed squares indicate Ohmic contacts to the 2DEG. (b- 
d) Characterization measurements: finite-bias stability plot 
of the QD (b), zero-bias conductance G of QPCl versus 
gate voltage Vqpci (c), and current- voltage characteristics 
of QPCl for different values of the gate voltage Vqpci close 
to pinchoff (d). Go is the conductance quantum. 



of the device. Aluminum surface Schottky gates were 
patterned on the GaAs/AlGaAs heterostructure by 
electron beam lithography, thermal evaporation and 
liftoff. The 2DEG is characterized by a density rig = 
2.26 X 10" cm-2 g^^^j ^^ mobility fi = 3.31 x lO^cm^/Vs. 
When a negative voltage is applied to surface gates, a 
central domain of area Au = 16 /im^ is defined by lat- 
eral confinement together with a QD and three QPCs. 
The measurement configuration is schematically shown 
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FIG. 2. (color online) . (a) Energy-level diagram for transport 
through a temperature-biased QD (all quantities are defined 
in the text), (b) Calculated zero-bias conductance G{e) for 
different temperatures of the leads: Ti = T2 = 250 mK (dot- 
ted line), Ti = T2 = IK (dashed line), Ti = 250 mK and 
Ta = iK (full line). 



in the same Fig. 1(a). Experiments were carried out in a 
■^He cryostat down to 250 mK. The domain was heated by 
driving a dc current Ih with a floating source through the 
QPCl-domain-QPC2 circuit while the zero-bias differen- 
tial conductance was simultaneously measured across the 
series QPC3-domain-QD [T3]. 

The single-particle energy-level spacing associ- 
ated with the lateral dimensions of the domain 
[Tr'^h'^/{2m*AD) « 400 neV] is much smaller than the 
thermal energy {ksT ~ 21iieV at 250 mK, where fc^ is 
the Boltzmann constant and m* is the effective mass), so 
that electrons in the domain can be treated as a Fermi 
gas. The charging energy of the domain was also found 
to have negligible impact on the low-bias measurements 
at the temperatures here of interest. Figures l(b-d) 
show transport measurements performed on the QD and 
QPCl with all other gates set to ground. The QD geom- 
etry was chosen in order to maximize its charging energy 
Ec and its single-particle energy spacing 5E. In the 
diamonds with lowest occupation number we measured 
Ec = 1.5 meV and SE « 0.4 meV [see Fig. 1(b)]. All 
QPCs show clear zero-bias conductance steps, as shown 
in Fig. 1(c) for the case of QPCl. Fig. 1(d) shows QPCl 
current- voltage characteristics at different gate voltages: 
as the gate voltage approaches pinchoff curves become 
increasingly non-linear in the given current interval. 

Based on the framework of Ref. [10 , we developed a 
model describing the CB-peak lineshape when a finite 
temperature bias is applied across the QD. We assumed 



for the QD a spin-degenerate single-particle energy spec- 
trum and an orbital level spacing exceeding thermal en- 
ergy. As a consequence, only one orbital level contributes 
to conduction for a given peak and we can model the QD 
as an Anderson impurity [141 115] with energy level e and 
on-site interaction Ec- The dot is tunnel-coupled to ideal 
leads at chemical potential /ii (/i2) and temperature Ti 
(T2). Coupling is assumed to be weak, so that transport 
is dominated by sequential tunneling and finite-lifetime 
broadening effects can be neglected. The energy diagram 
of the modeled system is sketched in Fig. 2(a). In the 
limit of large charging energy {Ec ^ ^5^1.2) and close 
to the resonance /ii = /i2 = e, the current is given by 
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where fi{e) = [l + e 



{>i-tl,}/kBT, 



1 is the Fermi-Dirac 



distribution function of the i"^ lead, Ti is the tunneling 
rate between the dot and the i**^ lead and e is the elec- 
tron charge [IS]. Fig. 2(b) shows the calculated G versus 
e in the degenerate case for some representative values of 
Ti and T2, assuming Ec = 1-5 meV. When Ti = T2 tem- 
perature variations rescale peaks without affecting their 
shape (black, green lines). By contrast, when Ti ^ T2 
the line-shape is distorted: peak tails are primarily de- 
termined by the higher temperature while the peak body 
by the lower one (orange line). 

The presence of this qualitatively different behav- 
ior is very useful for a comparison with experimental 
results. In doing so, we shall assume that the fast 
electron-electron scattering rate [17] does establish a 
quasi-equilibrium regime for the confined electron sys- 
tem, so that the microdomain can play the role of the 
hotter lead. Figures 3(a,b) show the measured zero-bias 
conductance G versus gate voltage Vc at Ttath = 250 
mK. The different curves correspond to increasing heat- 
ing currents Ih driven through the QPC1-QPC2 circuit. 
Here, the device is tuned so that Rqpci,2,3 — 7AkQ, 
while the conductance of the QD is of the order of a few 
jiS. These experimental characteristics closely resemble 
the calculated curves of Fig. 2(b). In order to perform 
a quantitative comparison, the derivative of Eq. IT] was 
fitted to the experimental points. The QD barriers were 
tuned to be symmetric, therefore we put Fi = r2 in the 
least-squares fit. Conversion of the gate voltage Vc into 
the dot energy e was provided by the finite-bias stabil- 
ity plot [TT]. Finally, after normalizing the curve to the 
position and value of its maximum, the two lead tem- 
peratures are the only free parameters. Following this 
procedure we can associate two temperatures (Ti and 
T2 , with T2 > Ti ) to each measured conductance charac- 
teristic. Figure 3(c) shows the extracted values of these 
temperatures as a function of the injection current Ih. 
T2 (in red) steadily increases with Ih, while Ti (in blue) 
exhibits only a small change, its value remaining close 
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FIG. 3. (color online). (a,b) Conductance G versus gate volt- 
age Vg for increasing values of the heating current In'- 11-6 n A 
(squares), 20 nA (circles), 28.6 nA (triangles). The corre- 
sponding full lines are best-fits of the derivative of Eq. IT] to 
each dataset. The hotter and colder temperature thereby es- 
timated are: Ti = 426 mK,!!^ = 604 mK, Ti = 446mK,r2 = 
861 mK and Ti = 485mK,r2 = 1.15 K, respectively. The 
data are plotted on both linear (a) and semi-logarithmic scale 
(b); in (a) traces have been offset vertically for clarity, (b) 
Extracted temperatures Ti and T2 as a function of Ih- The 
hotter T2 pertains to the heated microdomain, the colder Tj 
to the external portion of the 2DEG. 



to the one measured in an open configuration where only 
the gates defining the QD are polarized. We hence assign 
T2{Ih) to the central domain and Ti{Ih) to the adjacent 
2DEG region. There is a ~ 150 mK difference between 
the cryostat base temperature and the electronic temper- 
ature in the absence of injection that we ascribe to power 
leaking through cryostat lines. The additional increase 
of Ti with Ih stems from the global device heating pro- 
duced by large injection currents. The present procedure 
allowed us to measure local temperature differences from 
about 150 niK up to about 1 K. When T2-T1 < 150 mK, 
the fitting algorithm fails due to the proximity of two 
symmetrical local minima in the parameter space. 

After demonstrating the functionality of this thermom- 
etry scheme, we turn to the investigation of heat relax- 
ation mechanisms in the microdomain. At steady state, 
the temperature T2 is the result of the following thermal- 
balance equation: 



E^o 



PCi + QqD + Qe-ph — 
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FIG. 4. (color online) (a) Scheme of the main contributions to 
the steady-state power balance, (b-c) Temperature T2 versus 
injected power Ph for different configurations of the device: 
Rqpci,2,3 = 7.4 m (b), Rqpci,2 = 20Mfi, Rqpcs = 1 Mfi 
(c). For each case, the full lines are a fit of the expressions 
given in the text, (d) Joint bi-logarithmic plot of (b) and (c). 



the cryostat base temperature (Tph — T^ath — 250 mK) , 
while Qqpci and Qqd are the heat currents flowing 
through the i-th QPC and the QD, respectively |,18] . 
A schematic drawing of these contributions is shown in 
Fig. 4(a). We now briefly discuss each term. 

The power flow Qqpci through each QPC can be cal- 
culated within the Landauer-Biittiker formalism P^ [20] • 
If the transmission coefficient is energy-independent in 
the relevant range (which also results in a linear current- 
voltage characteristic), the result can be evaluated an- 
alytically. For the biased QPCs it takes the form 



Qqpci — Qqpc2 



VhIh 

4 



[T^ - Tl] , where 



where Qe-ph is the power exchanged between electrons 
and the phonon bath assumed to be at equilibrium at 



Vh is the total voltage drop developed across the heat- 
ing circuit, Lq — ir^kg/Se^ is the Lorenz number and 
we made use of the fact that Rqpci — Rqpc2- The 
first term in this expression is reminiscent of Joule heat- 
ing: of the total power VrIh provided by the current 
source, one half is delivered into the domain, Ph = 
IhVh/2. The second term is related to heat leakage 
from the domain to the (colder) adjacent 2DEG region 
through quasiparticle tunneling. It is proportional to 
the difference of the squared temperatures, in agree- 
ment with the Wiedemann- Franz law |21j . For the un- 
biased QPC, only the latter term contributes so that 
QqPC3 = {Lo/2Rqpc3) [T2 - Tf ] . 

Let us now turn to Qe-ph, the power transferred 
through electron-phonon coupling. At low tempera- 
tures, heat exchange between the 2DEG and the three- 



dimensional phonon bath mainly relies on the scatter- 
ing with low-energy acoustic phonons. In particular, 
for III-V semiconductor alloys, piezoelectric coupling is 
expected to dominate over deformation-potential inter- 
action [21] at subkelvin temperatures. The energy- loss 
rate for the screened piezoelectric electron-phonon in- 
teraction in the clean limit is given by the expression 

Qe-ph = SAd ( T2 — Tp^ 1 , where the coupling constant 
E depends on the electron density and on host-crystal 
material parameters [53]. For our GaAs/AlGaAs 2DEG 
we obtain E « SOfWiim-^R-^ [51]. 

In light of the foregoing discussion, the temperature 
data of Fig. 3(c) are re-plotted in Fig. 4(b) as a func- 
tion of the injected power Ph, evaluated from the mea- 
sured Vh and I^. The data are well described by the 
power law Ph = A + B{T^ - T^), where Ti = 400 mK 
and the constant A accounts for spurious heating. The 
best fit of this expression yields A = (—0.15 ± 0.02)pW 
and B = (5.0 ± 0.1)pWK-2. Notably, the value for B 
agrees with the 4.9 pWK^^ predicted for the heat leak- 
age through three QPCs of equal resistance 7.4 kfi. In 
this regime, heat flow through the QPCs dominates over 
electron-phonon heat exchange. By increasing the resis- 
tance of the QPCs, one can make the former contribution 
negligible. The data plotted in Fig. 4(c) were taken for a 
device configuration in which the QPCl and QPC2 zero- 
bias resistances were set to 20Mri, while that of QPCS 
to IMJl. The continuous line is a fit to the data of the 
power law Ph — C + -/^(Tf — TAJ. The good agreement 
with the experimental data demonstrates the transition 
to a regime where heat relaxation relies on coupling to 
the phonon bath. This is manifest in the joint log-log 
plot of Fig. 4(d), where the two datasets are charac- 
terized by different exponents as well as by a different 
order of magnitude for Ph- From the least-squares fit 
we extract C = (-5 ± l)fW and D = (170 ± 20)fWK-^ 
Since Ai, = IGjim^ this yields E = (ll±l)fWiim-2K-5. 
The obtained value is about three times smaller than the 
theoretical prediction and about two times smaller than 
what was previously measured for a much larger electron 
reservoir |25j . The reasons for such discrepancy are not 
completely clear. For our case, possible sources of un- 
certainty are related to the estimate of the effective area 
relevant for the electron-phonon interaction and of the 
average electron density inside the reservoir. 

In conclusion, we have demonstrated local electronic 
thermometry based on a lateral QD, shown that heat ex- 
change between the 2DEG and phonons depends on tem- 
perature according to a power law consistent with the 
screened electron-acoustic phonon piezoelectric interac- 
tion, and provided a measurement of the electron-phonon 
coupling constant in the subkelvin regime. Such results 
are relevant to the implementation of nanodevices which 
rely on the manipulation of thermal distributions, e.g. in 
electronic refrigeration schemes based on QDs [S] UHl HZ] ■ 



The method here reported can be readily applied to other 
material systems characterized by different electronic and 
phononic properties, e.g., InGaAs alloys with high In con- 
tent [25] or semiconductor nanowires [25] . 
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